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NOliENCLATURE 


A,  B constants,  Eq . (2) 

b nozzle-to-wall  separation  distance 

Cp  specific  heat  at  constant  pressure 

f,  fj,  functions,  Eqs.  (7)  and  (A2) , respectively 

g acceleration  of  gravity 

H plume  source-to-ceiling  distance 

h heat  transfer  coefficient,  Eq . (1) 

h characteristic  heat  transfer  coefficient,  Eq . (4) 

k kinematic  momentum  flux  of  jet 

Pr  Prandtl  number 

Q energy  release  rate  of  fire 

Qj  enthalpy  flux  of  jet 

Q*  dimensionless  Q,  Eq . (4) 

q rate  of  heat  transfer,  Eq . (1) 

qc  rate  of  in-plate  radial  heat  transfer 

q^  rate  of  heat  transfer  to  insulation 

Reg  Reynold’s  number,  Eq . (4) 

Re^  Reynold’s  number  defined  in  [3,  5] 

r radial  distance  from  impingement  point 

*1/2  r jet  where  velocity  is  half  of  axial  value 

Tag  surface  temperature  for  adiabatic  ceiling 

Tamb  ambient  temperature 

Tg  ceiling  surface  temperature 

t time 

V velocity  profile  of  wall/ceiling  jet 

V „ maximum  of  V 

max 


-v- 


function,  Eq . (2) 

b divided  by  diameter  of  jet  orifice 
distance  from  wall 
constant,  Eq.  (A4) 

(largest)  z where  V/Vmax  = 1/2 
plate  emissivity 

rate  of  heat  transfer  retained  in  plate 
dimensionless  Ta(j,  Eq . (5) 
fraction  of  Q lost  by  radiation 
kinematic  viscosity 
Stephan-Bolzmann  constant 


THE  BUOYANT  PLUME-DRIVEN  ADIABATIC  CEILING  TEMPERATURE  REVISITED 


Leonard  Y.  Cooper  and  Anne  Woodhouse 
Abstract 

In  previous  works,  convective  heat  transfer  from  buoyant  plume-driven 
ceiling  jets  to  unconfined  ceilings  has  been  estimated  using  a formula  for  the 
temperature  distribution  below  an  adiabatic  ceiling,  Ta(j,  obtained  from  exper- 
imental data  in  the  range  0 <_  r/H  <0.7  (r  is  the  radial  distance  from  the 
plume,  H is  the  plume  source-to-ceiling  distance).  The  present  study  re- 
evaluates this  data,  and  develops  an  independent  estimate  for  Ta^.  The 
analysis  takes  account  of  the  effect  of  ceiling  surface  re-radiation,  and  use 
is  made  of  the  previously  established  similarity  between  plume/ceiling-  and 
jet/wall-driven  heat  transfer  phenomena.  The  latter  similarity  is  the  basis 
of  a correlation  of  recently  reported  free  jet-wall  jet  "recovery  temperature" 
data  into  a normalized  Ta(j  distribution.  All  of  the  analysis  leads  to  new 
formulae  for  estimating  the  convective  heat  transfer  to  ceilings  during  enclo- 
sure fires.  These  new  results  confirm  previous  formulae,  and  extend  them  into 
the  larger  range  0 _<  r/H  _<  2.2. 

1 . INTRODUCTION 


The  convective  heat  transfer  to  ceilings  during  enclosure  fires  can  be 
related  to  the  heat  transfer  to  unconfined  ceiling  surfaces  from  buoyant 
plume-driven  ceiling  jets  [1].  Results  of  unconfined  ceiling  experiments  [2) 
suggest  a relatively  simple  means  of  estimating  this  unconfined  ceiling  heat 
transfer  rate,  q,  to  an  unconfined  non-adiabatic  ceiling,  as 
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<1  - h (Tad  - V 


(1) 


Figure  1 depicts  the  phenomena  in  question  near  the  fire  plume.  In  Eq.  (1), 

Tad  is  defined  as  the  temperature  distribution  (a  function  of  radial  distance, 

r,  from  the  plume-ceiling  impingement  point)  established  at  the  surface  of  an 

adiabatic  ceiling  (i.e.,  perfectly  insulated  and  with  zero  absorbtivity  and 

emissivity)  by  the  gas  flow  from  the  plume  of  a fire  of  strength,  Q,  and 

effective,  fire-to-ceiling  distance,  H,  T is  the  instantaneous  lower  surface 

ceiling  temperature  distribution,  and  h is  the  heat  transfer  coefficient  based 

on  the  Ta(j  - Tg  temperature  difference.  In  the  experiments  [2]  the  transient 

thermal  response  of  unconfined  0.00159  m thick,  cold-rolled  steel  plate 

ceilings,  insulated  on  the  top  side  and  initially  at  ambient  temperature  were 

measured  when  heated  by  weakly  radiating  premixed  propane  fires  of  constant 

strength  between  1.17-1.53  kW  (H's  were  between  0.58  and  0.81  m) . Analysis  of 

the  data  which  led  to  the  Eq.  (1)  formulation  used  long-time,  "near"-ceiling 

gas  temperatures  to  estimate  T (t  + »),  which  was  then  used  as  a surrogate  for 

s 

the  Ta(j  distribution. 

The  original  h results  [2]  did  not  distinguish  between  variations  in  the 
characteristic  Reynold's  number,  Re^,  of  the  plume  (which  only  varied  in  the 
narrow  range  1.5(1 0^ ) < Re^  < 2.0(10^)  in  the  laboratory  experiments).  In  a 
subsequent  study  [3],  it  was  noted  that  inertial  forces  compared  to  buoyancy 
forces  are  generally  large  (i.e.,  characteristic  Richardson  numbers  are  small) 
in  plume-driven  ceiling  jets  for  the  most  interesting  range  of  small  to 
moderate  values  of  r/H.  For  this  range  it  was  conjectured  that  the  near- 
ceiling flow  and  heat  transfer  characteristics  of  interest  could  be  directly 
related  to  the  analogous  characteristics  of  wall  or  ceiling  jets  driven  by 


heated  or  unheated  free  turbulent  jets  (e.g.,  the  same  configuration  of  Figure 
1,  but  with  a turbulent  jet  replacing  the  fire  and  its  buoyant  plume).  This 
led  to  a proposed  correlation  for  h which  does  depend  on  the  impingement  point 
Reynold's  number,  Re^,  namely 

I A Re"1/2  Pr"2/3  (1  + B r/H) , 0 < r/H  < 0.2; 

H. 

(2) 

C Re”0,3  Pr“2/3  w(r/H),  0.2  < r/H 

n 

where  A and  C are  numerical  constants,  B = B(Re^) , h/h  is  a continuous 
function  of  r/H  and 


lim  w ~ (r/H)  ^*2  (3) 

r/H  ->■  ® 


In  the  above 


£ _ n r 1/2  „l/2  *1/3 

h * pamb  Cp  8 H Q 


R.h  - g1/2  H3/2  Q*1/3/v 


(4) 


Q*  - (1  - X )Q/(p  C T g1/2  H5/2) 
r amb  p amb 


where  Xr  is  the  fraction  of  Q which  is  lost  by  radiation  from  the  fire's 
combustion  zone  and  from  the  plume  itself,  i.e.,  (1  - Xr)  Q is  the  effective 
portion  of  Q which  actually  drives  the  plume's  upward  convection.  Also, 

Pamb>  Cp>  ^amb»  v>  an^  ?r  are  the  density,  specific  heat  at  constant  pressur-  , 
temperature,  kinematic  viscosity  and  Prandtl  number  of  the  ambient  environ- 
ment. 
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The  above  correlation  was  based  on  measured  heat  transfer  (to  an 
elevated-temperature  wall)  [5]  and  flow  properties  [10]  of  unheated  free  jet- 
driven  wall  jets  (where  Ta(j  = Tam^)  , and  on  the  application  of  a criterion  of 
equivalence  for  free  jets  and  buoyant  plumes  at  their  respective  sections  of 
wall/ceiling  impingement  (Appendix  of  [3]). 

With  the  data  of  [2],  data  from  the  similar  "reduced-scale"  fire-plume- 
driven  heat  transfer  experiments  of  [4] , and  unheated,  free  jet-driven  wall 
jet  heat  transfer  measurements  of  [5],  it  was  shown  in  [3]  that  with  the 
proposed  h correlation  of  Eq.  (2),  the  Eq.  (1)  formulation  would  provide 
reasonable  q estimates  for  Re^  at  least  in  the  range  0.73(104)-27(104) . 

(Full  scale  hazardous  fires  with  fire-to-ceiling  distances  of  2-3  m and 
convected  energy  fluxes  of  200-2000  kW  would  have  Re^'s  in  the  range 
20(104)  < ReH  < 50(104)  [1]). 

In  [3]  the  h values  from  the  [2]  and  [4]  experiments  were  derived  from 
Eq.  (1)  with  the  use  of  measured  values  of  q,  Ts,  and  the  long-time,  "near"- 
ceiling  gas  temperature  measurements  of  [2]  which  were  assumed  to  be  represen- 
tative of  Ta(j.  The  h values  from  the  free-jet  experiments  of  [5]  were  derived 
in  a similar  way,  but,  as  mentioned  earlier,  the  values  of  Tfl(j  in  this  case 
required  no  special  measurements  as  they  are  identically  equal  to  Tam^. 

2.  A RE-EVALUATION  OF  EXISTING  DATA 

The  present  investigation  involves  a re-evaluation  of  the  data  of  [2] 
within  the  context  of  the  Eq.  (1)  formulation.  The  goal  is  to  obtain  new 
estimates  for  Ta<j  and  h which  would  1)  be  consistent  with  an  Eq.  (2)-type  of 
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ReH-dependent  h;  2)  would  only  use  the  initial  transient  data  (t  ■ 1,  2,  3,  5 
and  7 min)  of  [2],  and  not  the  25-30  min,  long-time,  near  surface  gas  tempera- 
tures which  were  previously  used  to  represent  the  "true"  Ta(j  [2];  3)  would 
include  the  previously  neglected  effect  of  exposed  ceiling  surface  re-radia- 
tion  to  far-field  surfaces;  and  which  would  4)  use  recently  reported  data  of 
free  jet-driven  wall  jet  heat  transfer  experiments  to  extend  the  r/H  range  of 
validity. 

With  regard  to  item  2)  above,  there  are  two  reasons  for  avoiding  the  use 
of  the  long-time  measurements  in  a new  analysis  of  the  data  of  [2].  First, 
there  is  indication  from  the  transient  data  that  in-plate  radial  conduction  is 
important  even  at  steady-state  for  r/H  at  least  up  to  0.25.  Also,  a non-zero 
plate  emissivity  could  lead  to  non-negligible  plate  radiation  to  the  far- 
field.  Under  either  of  these  circumstances  the  steady-state  plate  response 
would  differ  from  that  of  an  adiabatic  ceiling.  Second,  all  the  reported 
transient  data  are  from  thermocouples  mounted  directly  on  the  plate  itself, 
whereas  the  long-time  temperature  data  in  reference  [2]  were  measured  in  the 
gas  at  a distance  below  the  plate  surface  of  0.0016  m. 

Consistent  with  the  above  remarks,  another  analysis  of  the  tabulated  data 
of  [2]  was  carried  out.  The  analysis  involved  application  of  energy  conserva- 
tion of  elemental  volumes  of  the  plate,  at  the  seven  locations  of  plate- 
mounted  thermocouples  (placed  at  .07  m radial  intervals  from  the  point  of 
plume-ceiling  impingement),  and  at  different  times  during  the  plate  heating. 
New,  dimensionless  Ta^  data  were  generated  from  the  expression 
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(5) 
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- T 
ad  amb' 


ad 


Tamb  Q 


*27T 


(Ts  Tamb^ 

wn 


Tamb  Q 


hTambQ 


JJJ  + qc  + AE  + eo(T4 


- T , 
s amb 


)] 


where  o is  the  Stef an-Bolzmann  constant  and  where  q^ , qc  and  AE  are  the 
respective  instantaneous  rates  of  heat  transfer  per  unit  area  of  plate 
conducted  and  lost  to  the  rear  insulation,  conducted  radially  outward,  and 
absorbed  and  retained  in  the  plate  itself.  The  right-hand  term  of  Eq.  (5) 
represents  the  lower  plate  radiation  to  far-field  ambient  temperature  surfaces 
which  are  taken  to  be  at  the  ambient  temperature,  and  e is  the  effective  plate 
emissivity . 


A review  of  the  data  of  [2]  indicates  that  by  seven  minutes  into  the 
experiments  a value  of  e greater  than  a few  tenths  would  have  led  to  a signi- 
ficant contribution  (on  the  order  of  35  percent  for  e = 1)  to  the  net  value  of 
the  bracketed  term  of  Eq.  (5).  It  would  therefore  appear  that  an  Eq.  (5)-type 
of  analysis  of  the  data  requires  an  estimate  for  the  plate's  e.  Unfortunate- 
ly, the  literature  has  little  to  offer  for  e values  of  cold-rolled  steel 
plate.  Indeed,  while  tables  of  the  e's  of  solid  surfaces  typically  provide 
values  of  0. 6-0.8  for  rolled  steel  plate  of  various  descriptions,  the  authors 
found  only  one  primary  reference  [6]  for  the  e of  cold-rolled  steel.  Such 
an  e is  reported  to  have  been  measured  at  the  remarkably  low  values  of  0.075 
(93°C)  and  0.085  (260°C).  (Reference  [7]  reports  these  results  of  [6] 
incorrectly  as  0.75-0.85.)  To  confirm  this  the  authors  had  the  e of  available 
samples  of  "old  and  uncleaned"  cold-rolled  steel  measured  at  the  National 
Bureau  of  Standards.  A room  temperature  value  of  0.12  was  obtained  [8].  As 
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will  be  seen  below,  the  new  analysis  of  the  data  of  [2]  suggests  that  the 
"true"  e of  the  plate  of  [2]  was  indeed  in  the  range  of  these  low  values. 


Except  for  e (assumed  for  now  to  be  unknown)  and  h,  all  terms  on  the 
right  hand  side  of  Eq.  (5)  are  available  from  tabular  data  presented  in  [2]. 
To  carry  out  the  Eq.  (5)  analysis  the  following  h distribution  was  used 


8.82  Re^1/2  Pr“2/3  [l  - (5.0  - 0.284  Re°*2  )(r/H)],  0 < r/H  < 0.2 


h/h  = 


(6) 


0.283  Re 


-0.3  „_-2/3 

H 


-1.2  (r/H  - 0.0771) 


Pr  (r/H)  ’ (V/HTT279 r ’ °*2  * r/H 


This  expression,  which  is  in  the  form  of  the  correlation  of  Eqs.  (2)— (3) , 
was  developed  by  using  the  free  jet/buoyant  plume  equivalency  relations 
proposed  in  Eqs.  (A9)  and  (A10)  of  [3]  and  by  a curve  fitting  of  the  free-jet- 
driven  wall  jet  h data  of  [5]  (7.3(10^)  < Re^  < 27(10^)).  Plots  of  the  h data 
of  [5]  and  the  Eq.  (6)  correlation  are  presented  in  Figures  2 and  3 (residual 
standard  deviations,  rsd,  of  .87  and  .04,  respectively)  for  r/H  values  within 
(0  < r/H  < 0.2)  and  outside  the  stagnation  zone,  respectively.^ 


^Note  that  the  data  of  [5]  were  incorrectly  plotted  in  Figures  6-8  of  [3]. 
Instead  of  Re^*^  Pr^^  h/K  and  Re^'^  Pr^^  h/K  as  indicated  on  the 
ordinates  of  Figures  6-7  and  Figure  8,  respectively,  the  actual  plots  are  of 
Re^’^  Pr^/^  h/K  and  Re^*^  Pr^^  h/K,  where  Re^  is  defined  in  [3,  5]  and 
is  approximately  Re^  = 0.422  Re^.  All  of  the  data  are  correctly  plotted  lu-r.* 
in  Figures  2 and  3. 
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The  new  AT^  values  were  computed  from  Eqs.  (5)  and  (6)  and  the  data  of 
Tables  IV-IX  of  [2]  for  different  assumed  values  of  e ranging  from  0 to  1. 
Figures  4 and  5 are  plots  of  these  data  for  the  e = 0 and  e ■ 1 cases , 
respectively. 


3.  A FORMULA  FOR  AT*d(r/H)  AND  AN  INDEPENDENT  INDICATION  OF 

THE  "TRUE"  VALUE  OF  e 


The  correspondence  between  the  flow  and  heat  transfer  characteristics  of 

plume-driven  ceiling  jets  and  heated  or  unheated  free  turbulent  jet-driven 

wall  jets  was  mentioned  earlier.  Relative  to  this  correspondence  which  was 

conjectured  and  supported  in  [3],  and  which  was  introduced  into  the  present 

analysis  via  the  Eq.  (6)  description  for  h,  important  new  experimental  results 

have  been  reported  recently  in  [9].  Using  data  from  that  work,  which  studied 

heated  jet-driven  wall  flows,  the  development  of  the  following  expression  for 
* 

the  normalized  AT  ^ distribution  was  presented  in  the  Appendix 

ATad(r/H)  = f(r/R)  1-1 . 10(r/H)°*8  + 0.808(r/H) 1,6  (?) 

AT*  (r/H-0)  1-1. 10(r/H)0*8  + 2.20(r/H)1,6  + Q.690(r/H)2,4 

ad 

This  result  is  plotted  in  Figure  6 together  with  the  (equivalent)  data  from 
which  it  was  derived. 

A 

To  use  an  Eq.  (7)  correlation  for  the  new  AT^  values  it  is  necessary  to 

determine  the  "best"  value  for  the  e,  AT  ,(r/H=0)  pair.  Toward  this  end  the 

ad 

* 

following  procedure  was  carried  out:  1)  compute  new  AT^  values  for  a given 

e as  discussed  above;  2)  obtain  the  corresponding  ATa(j(r/H=0)  which  provides  a 
least  squares,  Eq.  (7)  fit  to  the  new  values. 
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Carrying  out  the  procedure  led  to  the  following  results:  for  t in  the 

range  0 < e < 1,  the  values  of  AT*^(r/H=0)  and  the  rsd  increase  monotonically 
with  e over  the  relatively  narrow  ranges  from  8.59  and  0.56  to  9.75  and  0.78, 
respectively.  Thus,  the  present  correlation  is  consistent  with  "very  low" 
values  of  plate  emissivity  for  which  the  smallest  rds’s  are  obtained. 

The  above  procedure  favors  e = 0 as  the  "best"  or  "true"  value.  However, 
in  view  of  the  above-mentioned  results  of  [6]  and  [8] , 

AT*  (r/H  = 0)  = 8.70  (8) 

ad 

which  corresponds  to  e = 0.1  and  a rsd  of  0.57  would  appear  to  be  a preferred 

value  for  the  normalizing  factor  of  Eq.  (7).  Using  Eqs.  (7)  and  (8)  leads  to 

* 

a new  estimate  for  AT  ^ which  is  plotted  in  Figure  7 along  with  the  new, 

* 

e = 0.1,  AT  , data, 
ad 

Also  plotted  in  Figure  7 is  a correlation  from  [3]  of  the  measured,  late- 
time, near-ceiling  gas  temperatures  of  [2].  As  can  be  seen,  except  for  small 

* 

r/H  values  the  new  estimate  for  AT^  is  close  to,  and  slightly  above  the 

measurements.  Near  the  stagnation  point,  however,  the  measured  temperatures 

generally  suggest  a plate  temperature  distribution  with  relatively  large 

radial  gradients  and  with  an  r/H  = 0 value  which  actually  exceeds  the  new 
* 

AT^.  The  large  gradients  can  be  explained  by  the  non-adiabatic,  in-plate, 
radial  conduction  effect  discussed  earlier.  However,  near-ceiling  gas  temper- 
atures would  never  have  exceeded  Ta(j,  and  the  fact  that  they  appear  to  do  so 

* 

would  seem  to  disqualify  the  Eq.  (7)  estimate  of  AT^  near  the  stagnation 
point. 
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In  [2]  it  was  noted  that  "appreciable  error"  in  the  estimates  of  in-plate 

radial  conduction  (q^  of  Eq.  (5))  "has  to  be  expected"  near  r/H  - 0.  This 

fact  together  with  the  above  discussion  and  previous  remarks  suggest  1)  that 

the  long-time  near-surface  temperature  measurements  of  [2]  continue  to  be  the 
* 

basis  of  the  AT^  estimate  within  the  stagnation  zone,  albeit  with  some 

uneasiness;  2)  the  latter  estimate  be  continued  at  r/H  = 0.2  by  a new  Eq.  (7)- 
* 

type  AT  ^ estimate  which  would  be  used  outside  the  stagnation  zone;  and 

3)  that  the  selection  of  an  optimum  normalizing  factor  to  use  in  Eq.  (7)  be 

* 

made  without  consideration  of  those  new  AT&^  values  which  correspond  to  the 
two  smallest  radial  positions,  r = 0 and  r = 0.07  m,  of  the  experiments. 


The  above  ideas  were  implemented,  and  they  resulted  in  the  following 

* 

newly  recommended  estimate  for  AT  , : 

ad 


10.22  - 14.9  r/H,  0 < r/H  < 0.2; 
8.39  f(r/H),  0.2  < r/H 


(9) 


where  f(r/H)  is  defined  in  Eq.  (7). 

4.  CONCLUSIONS 

* 

The  new  AT^  estimate  of  Eq.  (9)  provides  an  experimentally  based 
extension  of  earlier  results  from  a maximum  r/H  value  of  0.7  to  a r/H  value  of 
2.2.  The  results  of  the  present  analysis  provide  independent  confirmation 
that,  except  for  some  uncertainty  in  the  stagnation  zone,  the  measured,  long- 
time, near-ceiling  gas  temperatures  of  [2]  accurately  represent  values  of  Ta(j, 
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and  that  ceiling  re-radiation  in  those  experiments  was  negligible  because  of 
very  low  effective  values  of  surface  emissivity.  Also,  the  results  of  this 
paper  are  consistent  with  and  provide  additional  support  for  the  notion, 
introduced  in  [3],  that  up  to  moderate  values  of  r/H  there  is  an  equivalence 
between  the  flow  dynamics  and  heat  transfer  characteristics  of  free  jet-driven 
wall  flows  and  buoyant  plume-driven  ceiling  jets. 

A recommended  procedure  for  estimating  the  convective  heat  transfer  from 
a fire  plume  to  an  unconfined  ceiling,  from  the  point  of  plume  impingement  to 
moderate  r/H  was  developed  in  [3].  This  procedure  is  still  recommended. 
However,  the  new  results  for  Tad  and  h,  developed  here  and  presented  in  Eqs. 
(6),  (7)  and  (9)  should  replace  the  corresponding  Eq.  (11),  (12),  (17),  (18) 
and  (20)  estimates  of  [3].  The  new  results  should  also  be  used  to  modify  Eqs. 
(5),  (6),  (16)  and  (17)  of  [1]  where  the  unconfined  ceiling  equations  are  used 
to  estimate  heat  transfer  to  the  confined  ceilings  of  real  compartment  fire 
scenarios. 
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APPENDIX 


Ta(j  for  the  Wall  Jet  Away  from  the  Stagnation  Zone 


Consider  a wall  jet  driven  by  a heated  turbulent  jet  impinging  on  a plane 
adiabatic  surface.  Let  be  the  constant  enthalpy  flux  of  the  jet  relative 
to  Tam^.  Assume  that  jet  temperatures  near  the  impingement  point  are  close 
enough  to  Tam^  so  that,  for  the  purpose  of  applying  conservation  of  energy  to 
the  wall  jet,  the  density  field  and  Cp  can  be  approximated  by  their  ambient 
values.  Then 


T - T 


Q.  = 2irrp  ,C  V (T 
j amb  p max'' 


ad 


- T 


amb 


) / 

o 


amb 


ad 


- T 


dz 


(Al) 


amb , 


max 


where  the  integration  is  over  the  thickness  of  the  wall  jet  at  the  radial 
position  of  interest.  In  the  above,  V is  the  distribution  of  the  radial 
component  of  velocity  and  Vmav.  is  its  maximum  value.  Ta(^  is  the  temperature 
distribution  at  the  wall  surface. 


Let  6 be  the  distance  from  the  wall,  at  the  outer  portion  of  the  wall 
jet,  where  V/VTnav  = 1/2.  Then  the  experimental  study  of  [10]  for  unheated 
jets  indicates  that  V/Vmax  profiles  at  different  r stations  are  similar  in  the 
sense  that  they  can  be  closely  approximated  by 


V/V  = f.(z/6)  for  r/b  > 0.75 
max  1 


(A2) 
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where  b is  the  effective  jet  nozzle-to-wall  separation  distance.  We  assume 
that  this  latter  result  is  applicable  even  for  heated  jets  since  momentum 
forces  will  typically  dominate  buoyancy  forces  in  the  zero  to  moderate  r/b 
values  of  interest.  Furthermore,  it  is  reasonable  to  assume  that  (T-Tamb^/ 
<Tad-T  amb^  can  a-*-so  t>e  approximated  as  a function  of  z/6  in  the  r/b  > 0.75 
range.  Using  the  latter  assumptions  in  Eq.  (1)  leads  to 

T - T = a Q./(p  ,C  V r6) , r/b  > 0.75 
ad  amb  j amb  p max 

] f-d(f) 

o \ ad  amb  / max 

Besides  Eq.  (A2),  other  results  in  [10]  of  interest  here  are 

vmax/(^/b)  = 1.32  (r/b)"1,1;  5/b  = 0.098  (r/b)0,9  (A5) 

where  Pam^K  is  the  (constant)  momentum  flux  along  axial  sections  of  the 
unheated  jet  prior  to  impingement. 

For  the  case  of  the  heated  jet  it  is  reasonable  to  expect  that  the 
unheated  jet  results  of  Eqs.  (A5)  would  still  be  valid  provided  K and  b are 
redefined.  This  could  be  done,  for  example,  by  drawing  an  equivalence  at  the 
jet/wall  impingement  point  between  a heated  jet  of  interest  and  the  unheated 
jets  for  which  Eqs.  (A5)  are  valid.  This  notion  of  equivalence  was  used, 
e.g.,  in  [3].  There,  a buoyant  plume  (i.e.,  the  far  field  of  a limiting, 
small-K  or  low-momentum  jet)  was  taken  to  be  equivalent  to  an  unheated  jet  at 
respective  points  of  wall/ceiling  impingement  if  the  total  mass  and  momentum 
flux  of  the  plume  and  jet  were  identical.  Applying  these  criteria  for  equiva- 


(A3) 

(A4) 
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lence  led  to  the  results  in  [3]  that  the  flow  of  a buoyant  plume-driven 
ceiling  jet,  at  least  up  to  moderate  values  of  r/H,  would  simulate  the  flow  of 
an  unheated,  free  jet-driven  wall  jet  if  b and  K in  the  results  of  [10]  were 
replaced  by  point-heat-source-driven  plume  parameters  as  follows 


K1/2  = 0.636  g1/2  H3/2  Q*1/3;  b = 1.17  H 


(A6) 


Following  the  above,  it  is  assumed  that  the  idea  of  equivalence  at 
impingement  between  the  K,  b pair  of  an  unheated  jet  and  related  parameters  of 
a heated  jet  has  some  generality.  Accordingly,  Eqs.  (A5)  are  used  in  Eq.  (A3) 
with  the  result 

— . , 7.73  a(r/b)-0*8;  r/b  < 0.75  (A7) 

V(p-bV  /K) 

For  the  case  of  a buoyant  plume,  the  equivalence  relationships  of  Eq.  (A6) 
together  with  the  substitution  of  (l  - Xr)Q  for  Qj  are  used  in  Eq.  (A7)  to 
yield 


11.8  a(r/H)"0,8;  r/H  > 0.88 


* 

A Normalized  AT  , Distribution 
ad 


(A8) 


Ta(j  data  for  heated  jet-driven  wall  jets  have  been  recently  acquired  and 
plotted  in  Figure  9 of  [9].  This  plot  is  in  the  form 


(r)  - T 
(0)  - T 


amb 


amb 


(A9) 
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where  *1/2  is  defined  as  the  radial  position  of  the  impinging  jet,  at  the 
plane  of  impingement,  where  the  jet  velocity,  in  the  absence  of  the  wall, 
would  have  dropped  to  1/2  of  its  value  on  the  jet  axis.  Plotted  in  this  way, 
the  data  (r/r^2  ^ 20.6)  correlate  well  for  all  test  cases  where  the  ratio  of 
jet  orif ice-to-wall-distance  to  jet  orifice  diameter,  Z,  exceeds  1 (i.e.,  for 
jets  of  different  Qj  and  with  Z’s  of  5,  10,  and  15).  Interpreting  Eq.  (A7)  as 
an  asymptotic  result  for  large  r/rj^*  the  following  curve  fit  to  the  Z > 1 
data  of  [9]  has  been  obtained  with  a residual  standard  deviation  of  0.02 


In  the  experiments  of  [9],  buoyancy  forces  compared  to  inertial  forces  in 
the  impinging  downward-directed  jet  were  reported  to  be  small  enough  so  that 
within  jet  axis  distances  of  interest,  centerline  velocities  varied  on  the 
order  of  a few  percent  at  most  when  the  jet  was  turned  upward.  Accordingly  it 
is  reasonable  to  assume  that  the  velocity  distributions  (e.g. , the  values  of 
r- 1 / 2 ) of  the  jets  at  impingement  would  be  closely  approximated  by  the  corre- 
sponding velocity  distributions  of  similarly  configured  but  unheated  free 
jets.  Equation  (A9)  of  [3]  suggests  the  following  equivalence  between  the 
r^/2  °f  an  unheated  jet  and  the  H (distance  from  the  source)  of  a buoyant 
plume 


l-0.199(r/r1/2)°*8  + 0.0258(r/r1/2) 1,6 


l-Q.199(r/r1/2)0*8  + Q.0640(r/r1/2) K6  + 0.00443  0/r1/2)2*4 


(A10) 


r1/2  = 0.109  H 


(All) 
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Using  this  in  Eq.  (A10)  suggests  the  normalized  AT  , distribution  for  plume- 

ad 

driven  ceiling  jets  which  is  presented  in  Eq.  (7).  This  together  with  the 
data  of  [6]  (made  equivalent  with  the  use  of  Eq.  (All))  is  plotted  in  Figure 

6. 
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Figure  1.  Heat  transfer  to  an  unconfined  ceiling 


10.00 
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Figure  2.  Plot  of  the  h data  of  15),  7.3(1(T)  < ReH  < 27(104),  and  the  estimate 
of  Eq.  (6)  - In  the  stagnation  zone,  Q < r/H  < 0.2 


20 


Figure  3.  Plot  of  the  h data  of  [5],  7.3(10  ) < Re^  < 27(10  ), 
and  the  estimate  of  Eq . (6)  - outside  the  stagnation 
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Figure  4.  Plots  of  the  new  AT  , data  for  e “ 0 (test 
test  2:  A,  test  3:^C  12] ) 
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Figure  5.  Plots  of  the  new  AT  , data  for  e 1 (test 
test  2:  A,  test  3:aaX  [2]) 
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Figure  6.  Plots  of  the  normalized  AT^  data  of  [9]  and  the 
estimate  of  Eq . (7) 
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Figure  7.  Plot  of  the  new  AT  , data  for  e = 0^1  (test  1:  0, 

test  2:  A,  test  3f  X [2]),  the  AT  estimate  of 

Eqs.  (7)  and  (8)  ( ),  and  the  correlation 

(Eqs.  (11)  and  (12)  of  13])  for  the  measured,  long-time, 
near-ceiling  gas  temperatures  of  [2]  ( ) 
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